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Abstract 

Results from the Suzaku X-ray broad-band observations of clusters of galaxies are summarized. Aiming at understanding 
the physics of gas heating/particle acceleration and the cluster dynamical evolution, we search for non-thermal hard 
X-ray emission from merging clusters, particularly A2163 and the Bullet Cluster, based on the Suzaku and XMM- 
Newtonj Chandra joint analyzes. The observed hard X-ray emission is well represented by the single- or multi-temperature 
thermal models including super-hot (kT ~ 20 keV) gas. On the other hand, no significant non-thermal hard X-ray 
emission has been detected. Together with the presently available literature, the hard X-ray properties have been 
studied for about 10 clusters with Suzaku. The present status on Suzaku measurements of non-thermal X-ray emission 
and the cluster magnetic field are summarized and compared with those from the RXTE, BeppoSAX, and Swift satellites. 
The future prospects are briefly mentioned. 
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1 Introduction 



According to the standard scenario of the structure 
formation in the Universe, clusters form via collisions 
and mergers of smaller groups and clusters. A cluster 
merger has a kinetic energy of the order of 10 65 erg. 
This is the most energetic event in the Universe since 
the Big Bang. If two such objects collide with each 
other, a huge amount of energy may be released and 
a certain fraction is expected to heat the gas and ac- 
celerate particles through shock waves, and induce 
bulk and turbulent gas motions. 

Signatures of cluster merging can be recognized in 
many ways. In the X-ray band, irregular morphology 
and complex temperature structure of gas show that 
the system is disturbed due to the past mergers. At 
the radio wavelengths, synchrotron emission extend- 
ing over a Mpc scale have been discovered from more 
than 30 clusters [12]. The existence of radio halo 
emission suggests that relativistic electrons are being 
accelerated in the intracluster space. Interestingly, 
there is a correlation between the non-thermal, radio 
synchrotron power (Pi. 4) and thermal X-ray luminos- 
ity (Lx) for merging clusters while relaxed clusters 
without radio halo lie in a region well separated from 
the merging clusters on the Pi. 4 — Lx plane [4]. It 
is suggested that generation of high-energy particles 
is connected to the dynamical evolution of clusters 
[5]. From the observational point of view, however, 
the direct link between the radio synchrotron halo 
and the non-thermal hard X-ray property is yet to 
be clarified. 

In X-rays, the same population of high-energy 



electrons are thought to interact with 3K CMB 
photons and then generate non-thermal Inverse- 
Compton (IC) emission. The IC emission in excess 
of the thermal emission is then predicted to be seen 
in the hard X-ray band (typically, above 10 keV) 
where the thermal emission diminishes. In addition, 
from the radio observation alone, we cannot separate 
the energy of magnetic fields and the energy of high- 
energy electrons. However, by comparing the radio 
and hard X-ray fluxes, the cluster magnetic filed is 
also estimated under the assumption that same popu- 
lation of relativistic electrons scatter off of CMB pho- 
tons [H]. Models for non-thermal bremsstrahlung 
emission caused by suprathermal electrons with en- 
ergies of 10-200 keV are presented in [31] . 

The existence of non-thermal IC hard X-ray emis- 
sion has been reported in nearby clusters. The 
archetype is the Coma cluster: the non-thermal hard 
X-ray flux has been measured by RXTE [29] and Bep- 
poSAX [TO]. Recent reports based on the broad-band 
X-ray observations with Suzaku 35 = and Swift [36 
did not detect any significant non-thermal compo- 
nent and the hard X-ray flux is reproduced by ther- 
mal models. The mismatch between several satellites 
is suggested to be reconciled if different sizes of field- 
of- views are took into consideration[TT] . 

Difficulty in firm detection of non-thermal hard 
X-ray emission from clusters seems to indicate that 
the IC emission is in fact very faint and/or that its 
spatial distribution make it difficult to be detected by 
the present instruments with limited imaging capa- 
bility. Another issue is spectral modeling of ther- 
mal emission component. At energies > 10 keV, 
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the X-ray emission from hot clusters is still domi- 
nated by thermal bremsstrahlung emission [351 126) . 
Thus detailed characterization of thermal emission 
spectrum including the multi-temperature compo- 
nents is indispensable to separate the non-thermal 
emission from the thermal one. As is discovered in 
RX J1347. 5-1145 [2B], the hard X-ray emission in the 
Suzaku/H&rd X-ray Detector (HXD) band is origi- 
nated predominately by the super-hot (kT ~ 25 keV) 
gas in the cluster. Such violent mergers often produce 
shock-heated gas with high (3> 10 keV) temperature, 
which need to be properly took into account in the 
thermal modeling in order to study high-energy pop- 
ulation in the intracluster space. 

2 Objectives 

We search for non-thermal hard X-ray emission from 
merging clusters using the Suzaku broad-band X-ray 
spectroscopy to reveal the origin of hard X-ray emis- 
sion and obtain new insight into the gas physics and 
the cluster dynamical evolution. We also estimate the 
cluster magnetic fields in clusters by comparing the 
hard X-ray and radio fluxes. Finally we compare the 
Suzaku results with those from other satellites and 
summarize the present status of hard X-ray study of 
clusters and briefly discuss the future prospect. 

We adopt a cosmological model with £7m = 
0.27, 57a = 0.73, and the Hubble constant Hq = 
70 kms" 1 Mpc -1 throughout the paper. Quoted er- 
rors indicate the 90% confidence intervals unless oth- 
erwise specified. 

3 Sample and the Suzaku data 

The 5th Japanese satellite, Suzaku [23], is equipped 
with the X-ray CCD imaging spectrometers (XIS; 
[20] ) and HXD 33] and enables the sensitive broad- 
band X-ray observations thanks to the low and stable 
background levels. Particularly, the PIN detectors of 
the HXD instrument are useful for the present study 
because they have achieved the lowest background 
level than other missions in the 10-60 keV range [9]. 
For the purpose of constraining the hard X-ray prop- 
erties of merging clusters, we have conducted Suzaku 
observations of several targets. We focus on X-ray 
bright sources having radio synchrotron halos, par- 
ticularly A2163 at z = 0.203, and the Bullet Cluster 
at z = 0.296 since they are considered to have under- 
gone a recent, violent merger. The Suzaku's exposure 
times are 113 ksec and 41 ksec for A2163 south and 
north-east regions, and 80 ksec for the Bullet Cluster. 

Figure 1 shows the XMM-Newton image of A2163 
with overlaid the SuzakufHXD field of views for 
the two pointing observations. A2163 is the bright- 
est Abell cluster hosting one of the brightest ra- 



dio halo[T3l [7]- The previous X-ray observations 
showed the presence of high temperature region in 
the north-east, which can be attributed to the merger 
shock J3[ [21]. In the hard X-ray energies, [30] re- 
ported the detection of non-thermal IC emission with 
a flux level of F NT = 1-1±J;£ x lO" 11 ergs^cur 2 
in the 20-80 keV band from the long RXTE obser- 
vations, while [Hj derived an upper limit of _Fnt < 
5.6 x 10~ 12 ergs" 1 cm" 2 from the BeppoSAX/PDS 
observations. Thus the Suzaku observation will fur- 
ther provide independent measurements of hard X- 
ray properties of the sample. 




Figure 1: XMM-Newton image of A2163 (grayscale) 
with overlaid the Suzaku HXD field of views (34' x 34' 
in FWHM) for the two pointing observations (the 
solid boxes). The XMM-Newton integration region 
for the global spectrum is indicated with the dashed 
circle. 

The Bullet Cluster exhibits extended radio halo 
emission [16] and the prominent strong shock feature 
[2"2"] . The hottest gas with kT > 20 keV exists in the 
region of the radio halo enhancement. Possible detec- 
tion of non-thermal hard X-ray emission has been re- 
ported; the non-thermal flux in the 20-100 keV band 
is (0.5 ±0.2) x 10~ n ergs" 1 cm" 2 by RXTE [28] and 
3At\i x 10~ 12 ergs" 1 cm" 2 by Swift [2J. Thus it is 
worth examining the existence of IC emission with a 
detailed Suzaku analysis. 

4 Analysis strategy 

To accurately measure the non-thermal X-ray emis- 
sion from clusters at hard X-ray energies (> 10 keV), 
we need: 

1. careful assessment of background components, 
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2. detailed modeling of thermal emission. 

For 1., the Suzaku HXD background consists of the 
Cosmic X-ray Background and instrumental Non-X- 
ray Background (NXB) and is dominated by the lat- 
ter. We estimate the systematic error of NXB to 
be 2% by analyzing the HXD data during the Earth- 
occultation and quote the 2% systematic error (la) in 
the spectral analysis of A2163 and the Bullet Cluster. 
This is consistent with SUZAKU-MEMO 2008-03 by 
the instrument tearrfj]- 

For 2., since the merging clusters can have 
complex, multi-temperature structure including very 
hot thermal gas that emits hard X-rays, we apply 
single-, two-, multi-temperature thermal models to 
the Suzaku spectra. As detailed later, the multi- 
temperature model is constructed based on the anal- 
ysis of XMM-Newton or Chandra data. The Suzaku 
and XMM-Newton/ Chandra joint analysis allow us to 
take advantage of Suzaku's spectral sensitivity in the 
wide X-ray band and the XMM-Newton/ Chandra's 
high spatial resolution as demonstrated by [26] . 



5 Results 
5.1 A2163 



With Suzaku/RXD, 
from the cluster up 
the CXB and NXB 



we detected X-ray emission 
to 50 keV. With subtracted 
components, the 12-60 keV 
source flux is measured to be 1.52 ± 0.06 (±0.28) x 
10~ n ergs -1 cm~ 2 . Here the first and second errors 
indicate the 1-er statistical and systematic uncertain- 
ties. The l-o" systematic error of the flux is estimated 
by changing the normalization of the NXB model by 
±2%. Thus the detection of hard X-ray emission is 
significant at the > 50 level even if the systematic 
error of NXB is considered. 

Next we performed the joint XMM+ Suzaku/HXD 
spectral analyses under i) the single-temperature, ii) 
two-temperature, and iii) multi-temperature (multi- 
T) models. For i) and ii), the XMM-Newton PN and 
MOS spectra were accumulated from the r = 10' cir- 
cular region as shown in Figure 1. 

For i), the XMM+ Suzaku broad-band spectra in 
the 0.3-60 keV band can be fitted by the APEC ther- 
mal emission model: the resultant gas temperature 
and the metal abundance are kT — 13.5 ± 0.5 keV 
and Z — 0.29 ±0.10 solar, respectively, and x 2 /d.o.f. 
is 1249/1180. Thus the hard X-ray emission is likely 
to be dominated by the thermal emission. For ii), 
we have not find any significant improvement of the 
fits in comparison to the case i). On the other hand, 
there is obviously a complex temperature structure 
as shown by the previous X-ray observations, we at- 
tempted to construct the multi-T model for the ther- 



mal emission below. 

iii) The 10' spectral region of the XMM data is 
divided into 2' x 2' grids and the single-component 
APEC model is assumed for each grid. Fig- 
ure 2 shows the total model of multi-T components 
(namely the sum of the single APEC models for the 
grids) fit to the Suzaku/HXD data. We find that 
this multi-T model gives an acceptable fit to the 
HXD spectrum. Note that the model includes a 
kT ~ 18 keV APEC component for the north-east 
"shock " region, which is in agreement with the pre- 
vious XMM result [3] . The X-ray luminosity of this 
hot component is 5 x 10 44 ergs _1 , which contributes 
to the hard X-ray emission by about 15%. 




Energy (keV) 

Figure 2: Upper panel: the Suzaku/HXD spectrum of 
A2163 (the crosses) fitted with the multi-T + power- 
law model (the solid step function) . The contribution 
from the spectral components, i.e., multiple temper- 
ature components (the dotted lines) and the non- 
thermal power-law model (the dotted lines) are in- 
dicated in the figure. Bottom panel: the residual of 
the fitting is shown with the crosses. 

In the case of iii), the additional power-law com- 
ponent for non-thermal IC emission does not signif- 
icantly improves the fit to the observed HXD spec- 
trum. Assuming the multi-T+ power-law model with 
T = 2.18 (i.e., the same index in radio [7]), the 90% 
upper limit on the non-thermal emission is derived 
as 9.4 x 10~ 12 ergs -1 cm~ 2 . This gives a 3-times 
stronger constraints than the upper limit obtained 
from RXTE [SOj though the present result is con- 
sistent with RXTE and BeppoSAX 6 within their 
errors. 

More details of data reduction and discussions are 
to be described in a forthcoming paper (N. Ota et al. 
in preparation). 

5.2 The Bullet Cluster 

Following similar methods mentioned in H5.ll we an- 
alyzed the Suzaku and Chandra data of the Bul- 
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let Cluster. Our preliminary analysis showed that 
the global cluster spectrum taken with Suzaku is 
well represented by a two-component thermal model. 
On the other hand, the non-thermal emission is 
not significantly detected, yielding the 90% upper 
limit on a r = 1.5 power- law component to be 
- lCT 11 ergs^cuT 2 (20-100 keV). Note that the 
derived upper limit depends on the assumed photon 
index. 

The observed broad-band Suzaku spectrum in the 
1-50 keV band can be also fit by the multi-T model. 
Here the multi-T model is constructed by analyzing 
the deep Chandra observations. We did not confirm 
the previous IC detections with RXTE and Swift and 
the present results support the thermal origin of the 
hard X-ray emission. Under the multi-T model, the 
flux of hot (kT > 13 keV) thermal component is esti- 
mated to be Fhot ~2x 10 -12 ergs -1 cm -2 and com- 
parable to the non-thermal flux reported by [2] , while 
the IC flux is suggested to be even lower than .Fhot- 
More details of the analysis including assessment of 
systematic uncertainty due to the Suzaku- Chandra 
cross-calibration will be presented in Nagayoshi et 
al. in preparation. 



6 Discussion 

6.1 Non-thermal hard X-ray emission 

For the two hot clusters, A2163 and Bullet, the hard 
X-ray spectra taken with Suzaku is likely to be dom- 
inated by thermal emission, giving stronger limits 
on the non-thermal IC emission. We also note that 
the super-hot (~ 20 keV) gas exists in both clus- 
ters whose contribution is not negligible in the hard 
X-ray band. Thus this component should be prop- 
erly modeled in order to accurately measure the non- 
thermal component. Furthermore, such hot gas is 
over-pressured and thus thought to be short-lived, 
~ 0.5 Gyr [26]. Therefore the existence of super-hot 
gas supports the scenario of recent merger. 

Given the link between the radio flux and X-ray 
luminosity (see JQ), the generation of high-energy 
particles in intracluster space is likely to be connected 
with the merging events. Thus, is there any rela- 
tionship between Lx and IC flux (Tnt ) or the gas 
temperature (kT ) and i*NT ? 

With Suzaku, non-thermal hard X-ray emission 
has been constrained in about ten clusters. In Fig- 
ure 3, i*NT is plotted as a function of kT. The Suzaku 
results for nearby clusters as well as RX J1347.5- 
1145 were taken from the literature (T31 [T7| |T5l |55[ 
|33 [251 H [26]. The RXTE/ BeppoSAX/ Swift results 
[HI [U [36] are also quoted for comparison. Note 
that the non-thermal flux depends on modeling of 
thermal component as well as the assumption of 
the photon index for non-thermal component. For 



Suzaku, non-thermal emission has not been detected 
and the upper limit was derived for the nearby ob- 
jects including the Coma cluster, A2319 etc. Hence 
there is no clear relation on the Pnt — T plane. 

6.2 Cluster magnetic field 

Using the observed radio flux, S syn , and a relation 
Sic/S'syn = Ucmb/Ub [HI GZ] , we can infer the 
strength of the cluster magnetic field. In the case 
of A2163, adopting the radio flux S sya — 155 mJy at 
1.4 GHz 7 and the IC upper limit of Sic < 0.26 ,uJy 
at 12 keV from the this work, we obtain B > 0.09 fiG 
for the multi-T + power-law (r = 2.18) model. For 
the Bullet Cluster, the magnetic field is estimated to 
be B > 0.06 /iG under the two-temperature APEC 
+ power-law (r = 1.5) model. Thus the lower limit 
on the order of ~ 0.1 — 1/iG has been obtained in 
most of the clusters observed with Suzaku. 



-RXTE 

-BeppoSAX 

-Swift 



-Suzaku 




log(kT [keV]) 

Figure 3: Non-thermal IC hard X-ray flux in ten clus- 
ters measured by Suzaku. The results from RXTE, 
BeppoSAX, and Swift are also shown for comparison. 
See 36] for references. 



7 Summary 

In order to search for the non-thermal hard X-ray 
component and reveal the origin of high-energy emis- 
sion and the cluster dynamical evolution, the Suzaku 
broad-band X-ray observations of hot clusters hosting 
bright radio synchrotron halo, A2163 and the Bul- 
let cluster, have been performed. The Suzaku/HXD 
spectra for these two clusters are well represented by 
single- or two-temperature thermal models. Since the 
merging clusters can often have shock-heated, super- 
hot (kT ~ 20 keV) gas and its contribution is not 
negligible in the interpretation of the hard X-ray ori- 
gin, the determination of thermal component to high 
accuracy is indispensable. This motivated us to carry 
out the joint Suzaku+XMM-Newton/ Chandra spec- 
tral analysis under the multi-T model. As a result, 
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the observed Suzaku spectra are consistent with the 
multi-T models and we did not find any significant 
non-thermal IC emission from A2163 and the Bullet 
cluster. 

The properties of hard X-rays and the cluster 
magnetic fields have been studied for about ten clus- 
ters so far. There is no significant detection of IC 
emission by Suzaku, giving the upper limits on the 
level of ~ 10~ n ergs -1 cm~ 2 . Thus no clear rela- 
tionship between the non-thermal X-ray and radio 
properties or the IC luminosity {L\c)-T is seen. To 
further constrain the physics of gas heating and parti- 
cle acceleration in clusters, it is important to increase 
the number of cluster sample by applying the present 
analysis method to other clusters. 

Given relativistic particles are localized in small, 
shock regions, the present instruments with limited 
imaging capability may fail to identify the existence. 
If it is the case, the advent of hard X-ray imagers will 
benefit the cluster study. Now NuSTAR [14] is suc- 
cessfully in orbit and taking focused images of the 
high-energy X-ray sky. The ASTRO- H satellite is 
scheduled to be launched in 2014 [33] and will carry 
state-of-the-art instruments to realize sensitive X-ray 
observations in the 0.3-600 keV band. The Hard 
X-ray Imagers on ASTRO-H will have an effective 
area comparable to NuSTAR. These new hard X-ray 
imagers will enable more accurate high-temperature 
thermal component and identification of particle ac- 
celeration site to get higher signal-to-noise ratios. 
This will lead us to detect the IC emission to a level 
~ 2 orders of magnitudes lower than the present 
limit. In addition, the high-spectral resolution of 
X-ray micro-calorimeter will enable direct measure- 
ments of bulk/turbulent gas motions [24], Therefore 
ASTRO-H can measure the non-thermal energies in 
the form of kinetic gas motions (turbulence, bulk gas 
motion) and relativistic particles, leading us to draw 
a more comprehensive picture of the cluster structure 
and evolution. 
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